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The  analysis  of  the  stated  problems  involves  a  considerable  mathematical 
apparatus,  and  all  details  cannot  be  contained  in  a  relatively  short  fina] 
(report.  In  particular,  this  applies  to  the  analysis  of  the  beam  modes. 

On  the  other  hand,  a  fairly  complete  account  is  presented  of  the  backward 
wave  instability  problem. 

The  backward  wave  interaction  analysis  has  been  made  with  the  assumption 
of  confined-flow-focusing.  It  appears  possible  to  extend  the  analysis  to 
Brillouin-flow-focusing  with  uniform  or  periodic  magnetic  fields. 
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PREFACE 


This  final  technical  report  was  authorized  by  Dr.  T.  Wessel-Berg  and 
prepared  by  Electron  Dynamics  Division,  Hughes  Aircraft  Company, 
Torrance,  California,  on  Contract  F30602-79-C-0028  for  Rome  Air 
Development  Center,  Griff iss  Air  Force  Bat-i,  New  York.  It  summarizes 
the  results  of  a  continuation  effort  on  the  "Development  of  a  Large 
Signal  Computer  Theory  for  TWT,"  based  on  polarization  variables  that 
was  reported  on  Contract  F30602-77-C-0221.  This  continuation  effort 
describes  the  implementation  of  variable  circuit  parameters  along  the 
tube  axis,  and  provides  a  stability  analysis  of  backward  wave  inter¬ 
action,  based  on  the  assumption  of  confined-flow-focusing. 

This  effort  has  been  initiated  in  December  1978  and  was  completed  in 
March  1980.  Joseph  Polniaszek  was  RADC  Froject  Engineer. 


I 


t 


TABLE  OF  CONTEXTS 


Section 

1.0  INTRODUCTION 

2.0  THE  TREATMENT  OF  TAPERED  CIRCUIT  PARAMETERS  2 

2.1  Introduction  2 

2.2  Mathematical  Approach  2 

3.0  RELATIVISTIC  CORRECTIONS  4 

4.0  FUNDAMENTAL  AND  SPACE-HARMONIC  INTERACTIONS  IN 

THE  TUT  6 

4.1 

4.2 

4.3 

4.4 

4.5 

4.6 


4.6.1  Left-Handed  Helix  23 

4.6.2  Right-Hsnded  Helix  24 

4.7  Compact  Formulation  of  the  Full  Set  of  RF 

Equations  24 

4.7.1  Left-Handed  Helix  Section  25 

4.7.2  Right-Handed  Helix  Section  26 

4.8  Special  Cases  and  Approximations  27 

4.9  The  Transmission  Matrix  of  a  Uniform  Helix 

Section  29 

4.10  The  Over-all  Transmission  Matrix  of  the 

Multisection  TWT  32 

4.11  Two-Port  Terminal  Description  of  the 

Multisection  TUT  34 

5.0  BACKL’ARD  WAVE  INSTABILITIES  38 

5.1  The  Inherent  Backward  Wave  Instability  38 


Interactions  of  Beam  and  Circuit  in  the 

Backward  Modes  6 

The  Architecture  of  Forward  and  Backward 

Wave  Interactions  10 

The  Matching  Problem  12 

Dispersion  Relations  for  the  Forward  and 

Backward  Waves  16 

Specification  of  the  RF  Variables  in  the 

Backward  Traveling  Wave  20 

Summary  of  Normal  Mode  Amplitudes  23 


Bl 


v 


TrticsDiao  am  same-hot  filmd 


TABLE  OF  CONTENTS  (CONTENTED) 


Section  Page 


5.2  Resolution  in  Forward  and  Backward  Components 

in  the  Over-all  Tranmission  System  40 

5.3  Proper  Procedure  for  Determination  of  the 

Inherent  Backward  Wave  Instability  45 

5.4  The  Effect  of  Mismatched  Input  and  Output 

Terminals  on  Backward  Wave  Instabilities  48 

5.4.1  Forward  and  Backward  Characteristic 

Impedances  48 

5.4.2  The  Reflections  at  the  Input  and 

Output  Terminals  50 

6.0  PLASMA  REDl'CTION  FACTORS  AND  BEAM  COUPLING  COEFFICIENTS 

FOR  FORWARD  AND  BACKWARD  MODES  56 

6.1  Beam  Model  56 

6.2  Confined  Flow  57 

6.2.1  Plasma  Reduction  Factors  for  Confined 

Flow  59 

6.2.2  Coupling  Coefficients  for  Confined 

Flow  63 

6.3  Brillouin  Flow  65 

6.3.1  Plasma  Reduction  Factors  for  the 

Backward  Waves  (m  -  il)  in  Brillouin 
Flow 


67 


LIST  OF  ILLUSTRATIONS 


Figure 

4-1  Dispersion  diagram  of  left-handed  helix. 

4-2  Dispersion  diagram  of  right-handed  helix. 

4-3  Schematic  diagram  illustrating  the  circuit  and  beam 
harmonics  involved  in:  a)  forward  interaction, 
b)  backward  interaction  in  a  left-handed  helix. 

4-4  Schematic  diagram  of  forward  interaction,  and  the 
rf  variables  involved. 

4-5  Schematic  diagram  of  backward  interaction,  and  the 
rf  variables  involved:  a)  left-handed  helix, 
b)  right-har.ded  helix. 

4-6  Schematic  illustration  of  one  uniform  helix  section. 

4- 7  Schematic  configuration  of  a  multisection  TWT 

consisting  of  a  total  of  n  uniform,  but  different 
sections . 

5- 1  Illustration  of  instability  diagram  for  backward 

wave  oscillations, 

5-2  Sketch  illustrating  the  nature  of  backward  wave 
instabilities . 

5-3  Schematic  diagram  of  a  sectioned  TWT  with  arbitrary 

input  and  oi.cput  loads  Z,  ,  and  ZT„. 

LI  LN 

5-4  Sketch  showing  the  TWT  separated  into  two  regions, 
with  the  attenuator  load  Z  serving  as  output  load 
for  the  first  region  and  input  load  for  the  second 
region. 

5- 5  Schematic  illustration  of  the  condition  (5.36)  of 

unity  loop  gain  as  criterion  for  backward  wave 
oscillations . 

6- 1  Schematic  diagrams  of  the  assumed  distributions  of 

the  dynamic  variables  over  the  beam  cross  section. 

6-2  Plasma  frequency  reduction  factor  Rq  for  forward 
wave  interaction  and  confined  flow. 


8 


9 


11 

13 


14 

31 


33 

41 

47 

49 


49 


54 

58 

61 


vii 


LIST  OF  ILLUSTRATION'S  (CONTINUED) 


Figure  Page 

6-3  Plasma  frequency  reduction  factor  R+^  for  backward 

wave  interaction  and  confined  flow.-  62 

6-4  Beam-circuit  coupling  coefficient  f^  for  forward  wave 

interaction  (m  =  0)  and  confined  flow.  64 

6-5  Beam-circuit  coupling  coefficient  f+^  for  backward  wave 
interaction  (m  -  ±1)  and  confined  flow. 


66 


EVALUATION 


This  report  describes  a  novel  approach  to  large  signal  modeling 
of  a  helix  type  TWT.  It  embodies  the  benefits  of  increased  accuracy 
and  decreased  computation  time  for  initial  tube  modeling  and  develop¬ 
ment.  This  is  accomplished  through  the  use  of  a  transformation 
that  reduces  the  nonlinear  equations  describing  the  system  into  a 
set  of  linear  differential  equations  that  can  be  solved  using  classi¬ 
cal  techniques.  The  ultimate  benefit  of  this  program  will  be  de¬ 


creased  development  time  and  cost  through  a  reduction  in  the  com¬ 
putation  time  and  also  the  "cut  and  dry"  experimental  techniques 


used  presently. 
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1 . 0  INTRODUCTION 

This  report  presents  continuation  efforts  in  the  large  signal  computer 
theory  described  in  report  No.  W-07587.  It  represents  a  further  investi¬ 
gation  of  special  effects  and  conditions  of  importance  in  TWT  design.  In 
particular,  it  treats  two  main  areas:  Implementation  of  variable  circuit 
parameters  along  the  tube,  and  backward  wave  instabilities. 

The  first  of  these  main  areas  involves  the  handling  of  velocity  tapers, 
impedance  tapers,  and  variable  attenuator  sections,  all  in  any  prescribed 
fashion. 

The  second  area  is  concerned  with  analysis  of  backward  wave  instabilities 
under  the  same  conditions  of  variable  parameters  along  the  tube.  In 
particular,  the  theory  is  applicable  for  segmented  TWTs,  i.e.,  tubes  with 
several  sections  of  otherwise  constant  circuit  parameters.  Such  con¬ 
figurations  are  known  to  reduce  backward  wave  instabilities,  but  the 
underlying  theoretical  background  and  the  modeling  have  been  inadequate. 
In  particular,  the  symmetries  and  coupling  relations  of  beam  and  cir¬ 
cuit  modes  have  not  been  sufficiently  well  understood.  The  present 
work  is  an  effort  to  contribute  to  a  better  understanding  of  these 
effects  and  provide  improved  models  and  mathematical  procedures  for  use 
in  numerical  analysis  of  backward  wave  instabilities. 

One  part  of  the  investigation  is  the  analysis  of  the  properties  of  the 
beam  modes  involved  in  backward  wave  interactions.  These  have  different 
symmetries  from  the  usual  circularly  symmetric  modes  involved  in  for¬ 
ward  interaction,  and  need  special  consideration. 

The  analysis  of  the  stated  problems  involves  a  considerable  mathematical 
apparatus,  and  all  details  can  not  be  contained  in  a  relatively  short 
final  report.  In  particular,  this  applies  to  the  analysis  of  the  beam 
modes.  On  the  other  hand,  a  fairly  complete  account  is  presented  of 
the  backward  wave  instability  problem. 
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2.0  THE  TREATMENT  OF  TAPERED  CIRCUIT  PARAMETERS 


2.1  INTRODUCTION 

If  the  circuit  parameters  are  varying  along  the  z-axis,  the  model  and  the 
basic  formulation  must  satisfy  certain  requirements  which  can  be  stated 
as  follows:  (i)  the  basic  circuit  and  beam  rf-variables  must  be  chosen 
such  that  they  are  continuous  along  the  axis,  even  if  the  circuit  param¬ 
eters  are  discontinuous,  (ii)  the  beam-to^circuit  coupling  must  be 
appropriately  formulated  in  terms  of  wave  length- dependent  coupling 
coefficients  rather  than  frequency-dependent  coefficients. 

The  requirement  (i)  serves  to  facilitate  the  mathematical  formulation 
of  problems  having  z-dependence  in  general.  The  requirement  (ii)  is 
necessary  for  handling  cases  of  high  circuit  loss,  i.e.,  attenuators  in 
general.  In  such  cases  the  wavelength  beam-to-circuit  coupling  is 
specified  by  an  attenuated  and  increasingly  fast  propagation  factor 
entering  the  argument  of  the  coupling  coefficient,  which  then  itself 
becomes  complex.  This  implies  that  the  coupling  is  no  longer  pure 
capacitive  but  has  a  real  component  as  well. 

In  the  polarization  model  used  in  the  present  work  both  requirements  (i) 
and  (ii)  stated  above  are  fully  accounted  for. 

2.2  MATHEMATICAL  APPROACH 

The  following  is  a  formal  exposition  of  the  mathematical  procedure.  In 
the  first  part  of  this  work,  reported  earlier,  we  established  the  concept 
of  the  transmission  matrix  T  for  a  uniform  helix  section.  The  trans¬ 
mission  matrix  relates  the  rf  variables  at  the  input  and  output  ends  of 
the  uniform  section  in  the  following  way: 

a(Zp)  =  T (p)  a(Zp_1)  ,  (2.1) 


2 


where  a(z  )  is  a  column  vector  containing  the  rf  variables  of  all  the 
~  P 

frequency  components.  The  transmission  matrix  implicitly  contains  large 
signal  effects,  in  that  the  off-diagonal  submatrices  are  nonzero.  These 
represent  the  nonlinear  coupling  terms  between  the  various  harmonics. 

If  the  circuit  parameters  vary  over  a  certain  length  of  the  TWT,  this 
length  is  subdivided  into  a  sufficiently  large  number  of  subsections. 

In  each  of  these  the  parameters  are  considered  to  be  constant  and 
independent  of  z. 

f 

The  over-all  transmission  matrix  for  the  tapered  section  is  then  given 
by  the  product  of  the  separate  transmission  matrices  of  the  p 
subdivisions. 


T  =  T (p)  T(p  -  1)  .  .  .  T ( 2)  T(l)  (2.2) 

The  formula  applies  for  the  practical  computational  case  in  which  the 
number  p  is  finite.  The  exact  formula  is 

T  =  Urn  [T (p)  T(p  -  1)  .  .  .  T(2)  T(l> ]  (2.3) 

p  -V  00 

This  simply  means  that  the  subdivisions  have  to  be  made  sufficiently 
small  to  simulate  the  real  situation.  In  practice,  this  is  no  problem, 
since  the  computer  program  is  sufficiently  general  to  accommodate  the 
indicated  procedure.  It  is  more  a  question  of  organization  of  input  data 
for  the  computer. 

Tn  the  small  signal  domain,  the  over-all  transmission  matrix  is  diagonal 
on  a  submatrix  basis,  and  the  computations  are  correspondingly  simple,  in 
that  each  harmonic  components  can  be  treated  independently  of  the 
others  . 


3.0  RELATIVISTIC  CORRECTIONS 


The  relativistic  corrections  that  have  to  be  implemented  at  high  dc 
voltages  have  been  introduced.  These  corrections  are  nil  expressed  in 
terms  of  the  usual  relativistic  factor  y,  defined  by: 


The  relativistic  effects  appear  in  the  following  parts: 


(3.1) 


1.  The  relativistically  correct  dc  velocity  v^  is  obtained 
from  the  dc  voltage  by: 


Ii-V 
2  m  0 
o 


V 

l  +iJL_0 
2  m.  2 
0  c 


mn  2 
0  c 


(3.2) 


The  last  fraction  represents  the  relativistic  correction,  beinj 
unity  if  VQ  is  small. 


2*  The  longitudinal  relativistic  mass  m  is  increased  by  the 

3  e' 

factor  1 /y  . 


m 


£ 


=  m. 


1 

3 

Y 


(3.3) 
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3. 


4. 


The  transverse  relativistic  mass  m  is  increased  by  the 
factor  1/y. 


m 


t 


(3.4) 


The  plasma  frequency  oi 


P 


is  reduced  by  the  factor  y 


3/2 


01 

P 


(3.5) 


With  the  indicated  corrections  introduced  into  the  program 
at  the  appropriate  places  the  results  are  relativistically 
correct  for  small  signals  and  approximately  correct  for  large 
signals . 
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4.0  FUNDAMENTAL  AND  SPACE -HARMONIC  INTERACTIONS 

IN  THE  TWT 


The  main  objective  of  the  present  work  is  to  analyze  backward  wave 
instabilities.  But  in  order  to  accomplish  this  task  it  has  been  neces¬ 
sary  first  to  establish  a  sufficiently  detailed  physical  model  together 
with  a  sufficiently  simple  mathematical  formulation  of  the  over-all 
forward  and  backward  interactions  taking  place  in  a  multisection  TWT. 

f 

4.1  INTERACTIONS  OF  BEAM  AND  CIRCUIT  IN  THE  BACKWARD  MODES 

It  is  necessary  to  distinguish  between  left-  and  right-handed  helices 
because  of  possible  inclusions  of  pitch  reversals  in  the  helix  struc¬ 
ture.  The  rf  fields  of  left-  and  right-handed  helices  are  given  by  the 
fallowing  expansions  in  space-harmonic  components: 


Ur,  6,  z)  -  £  MO 


m 


m=-oo 


Ur,  0,  z)  -  Y  K (t)  c‘j6”Z  eJ"°, 

Immai  m 


m=-c° 


In  these  expansions  the  exponent  is  given  by: 


(4.1) 


(4.2) 


B 


m 


(4.3) 


where  d  is  the  helix  pitch. 
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The  dispersion  diagrams  of  the  left-  and  right-handed  helices  are  shown 
in  Figures  4-1  and  4-2,  respectively. 

The  figures  labeled  a)  and  b)  in  these  drawings  are  mirror  images  of 
each  other.  They  both  exist  because  the  helix  looks  the  same  in  either 
direction. 

From  a  study  of  the  diagrams  we  can  draw  the  following  conclusions 
relating  to  the  basic  interactions: 

1.  For  a  left-handed  helix  the  backward  wave  interaction 
involves  the  m  “  -1  space  harmonics. 

2.  For  a  right-handed  helix  the  backward  wave  interaction  involves 

*'  the  m  =  +1  space  harmonics. 

3.  The  fundamental  mode  interaction  is  not  affected  by  the  direc¬ 
tion  of  helix  pitch. 

4.  In  a  given  left-  or  right-handed  helix  section  only  one  of  the 
two  space  harmonics  m  =  -1  or  m  =  +1  is  involved.  This  is 
because  the  m  =  -l  and  m  =  +1  harmonics  are  always  in  opposite 
directions.  Therefore,  only  one  of  the  two  space-harmonic 
components  can  be  synchronized  with  the  beam,  namely  the  one 
with  positive  propagation  constant. 

The  same  considerations  do  not  apply  to  the  corresponding  space  har¬ 
monic  components  of  the  beam.  If  a  m  =  -1  mode  is  excited  in  the  beam 
through  backward  wave  interaction  in  a  left-handed  helix  section,  this 
particular  beam  symmetry  is  retained  through  a  pitch  reversal.  In  the 
reversed  pitch  section  the  m  =  -1  existing  in  the  beam  is  completely 
decoupled  from  the  circuit,  because  the  m  -  -1  circuit  mode  in  the 
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Figure  4-1  Dispersion  diagram  of  left-handed  helix. 

a)  Forward  wave  interaction  in  the 
m  =  0  mode 

b)  Backward  wave  interaction  in  the 
m  ■  -1  mode 
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Figure  4-2  Dispersion  diagram  of  right-handed  helix. 

a)  Forward  wave  interaction  in  the 
m  *  0  mode 

b)  Backward  wave  interaction  in  the 
m  *  +1  mode 


right-handed  helix  is  completely  out  of  synchronism  with  the  beam.  In’ 
such  a  pitch  reversal  the  effect  on  the  beam  is  that  it  suddenly  sees  a 
noninteracting  circuit,  all  referred  to  the  m  =  -1  beam  mode. 

Conversely,  if  we  started  out  with  a  right-handed  helix,  the  m  =  +1  mode 
is  the  appropriate  backward  mode.  A  pitch  reversal  reduces  the  beam- 
circuit  interaction  for  the  m  =  +1  mode  to  essentially  zero. 

The  points  brought  up  here  might  conceivably  be  of  importance  for  reduc¬ 
tion  of  backward  wave  instabilities  through  pitch  reversals.  Provided 
these  can  be  designed  with  negligible  fundamental  mode  reflections  over 
the  appropriate  frequency  band,  they  appear  to  represent  a  practical 
possibility.  As  far  as  backward  wave  interactions  is  concerned  the 
effective  length  of  the  helix  is  reduced  by  a  factor  of  two,  which  cer¬ 
tainly  would  serve  to  diminish  the  instability  problem. 

4.2  THE  ARCHITECTURE  OF  FORWARD  AND  BACKWARD  WAVE  INTERACTIONS 


Using  the  previous  discussion  as  a  guideline  let  us  establish  the  rele¬ 
vant  interactions  that  need  be  considered.  Figure  4-3  shows  a  schematic 
diagram  which  helps  sort  out  the  interactions  and  the  angular  sym¬ 
metries  of  the  modes  involved.  The  first  point  which  we  have  to  keep 
in  mind  is  that  the  circuit  alone  supports  two  independent  modes,  one 

traveling  in  the  forward  direction  with  group  velocity  v  >  0,  and  one 

b 

traveling  in  the  backward  direction  with  v  <0.  The  corresponding 

6 

interactions  with  the  beam  are  those  shown  in  the  lower  blocks  of 
Figure  4- 3a  and  4- 3b. 


The  second  point  which  we  must  keep  in  mind  is  that  all  circuit  har¬ 
monics  are  intimately  coupled,  so  that  the  excitation  of  one  particular 
harmonic  by  the  beam,  say  the  m  =  -1  harmonic  shown  in  Figure  4-3b, 
automatically  gives  rise  to  a  proport  innate  excitation  of  a  1 1  the  other 
circuit  harmonies. 


10 


G7M7 


CIRCUIT  FORWARD  MODE,  *g  >  0 


CIRCUIT 


Cl  RCUIT  BACKWARD  MODE,  vg  <  0 


CIRCUIT 


Figure  4-3  Schematic  diagram  illustrating  the  circuit  and 
beam  harmonics  involved  in:  a)  forward  inter¬ 
action,  b)  backward  interaction  in  a  left- 
handed  helix. 


The  third  point,  we  like  to  make  is  that  it  is  convenient  to  interpret  • 
the  circuit  voltage  and  circuit  current  as  the  fundamental  components, 
i.e.,  those  corresponding  to  the  m  =  0  mode,  both  in  forward  and  back¬ 
ward  interactions. 

The  relevant  interactions  and  rf  circuit  and  beam  variables  are  those 
summarized  in  the  schematic  diagrams  of  Figures  4-4  and  4-5.  The 
forward  interaction  shown  in  Figure  4-4  involves  only  the  circularly 
symmetric  modes  m  =  0  of  the  circuit  and  the  beam.  The  rf  variables 
are  the  normalized  circuit  voltage  v^,  the  circuit  current  i^,  the 
normalized  beam  velocity  Uq ,  and  the  normalized  displacement  Sq. 

This  mode  of  operation  is  the  regular  TWT  amplification  process,  which 
is  independent  of  the  direction  of  pitch. 

Thfe  backward  wave  interactions  depicted  in  Figure  4-5  involve  the 
m  =  ±1  beam  and  circuit  modes,  and  the  m  »  0  circuit  mode.  As  stated 
previously,  the  relations  between  the  m  =  0  and  m  =  ±1  circuit  space 
harmonics  are  fixed. 

4.3  THE  MATCHING  PROBLEM 


With  these  considerations  in  mind  let  us  consider  the  matching  problem 

at  some  discontinuity  along  the  helix,  such  as  a  change  in  pitch  angle. 

Which  rf  variables  are  continuous  across  the  discontinuity,  and  need 

be  considered?  Expressed  in  general  terms,  the  circuit  voltage,  circuit 

current,  beam  velocity,  and  beam  displacement  are  continuous.  Let  us 

first  consider  the  circuit  variables.  As  noted  previously  it  suffices 

to  express  these  in  terms  of  the  fundamental  (m  =  0)  component,  which 

consists  of  the  sum  of  the  forward  components  v^_  and  i^  (Figure  4-4) 

and  the  backward  traveling  components  v,  and  i.  (Figure  4-5).  The 

b  b 

latter  are  also  of  the  m  =  0  symmetry,  but  arc  caused  by  interactions 
of  the  m  =  -!  or  m  =  +1  components,  depending  on  the  sense  of  rotation 
of  the  helix,  as  already  shown  in  Figure  4-5. 


1.2 


Figure  4-4  Schematic  diagram  of  forward  interaction, 
and  the  rf  variables  involved. 


g;mi 


a)  b) 


Figure  4-5  Schematic  diagram  of  backward  interaction, 
and  the  rf  variables  involved: 

a)  left-handed  helix. 

b)  right-handed  helix. 
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The  rf  variables  of  the  beam  corresponding  to  different  m-values  are 
orthogonal,  and  presumably  are  not  coupled  by  a  change  in  helix  pitch. 
Accordingly,  it  is  necessary  to  match  the  various  rf  beam  variables, 
corresponding  to  different  m-values,  separately.  These  variables  are 
U_j  and  S  for  the  m  =  -l  symmetry,  U+j  and  S+j  for  the  m  +  1  sym¬ 
metry,  and  Uq  and  for  the  m  =  0  symmetry.  Altogether  we  have  eight 
matching  conditions  of  circuit  and  beam  variables.  If  pitch  reversals 
are  not  involved  at  any  position  along  the  helix,  these  are  reduced  to 
six  conditions. 

If  the  forward  mode  is  disregarded,  these  six  conditions  are  lurther 
reduced  to  four,  namely  continuity  of: 

Backward  circuit  voltage  v^ 

Backward  circuit  current  i, 

_  .  „  (4-4) 

Beam  velocity  U_j 

Beam  displacement  S  ^  , 
where  we  have  assumed  a  left-handed  helix. 

However,  with  the  forward  wave  disregarded,  there  are  only  three  inde¬ 
pendent  wave  solutions  available.  These  are  the  three  solutions  of  the 
third-order  backward  wave  dispersion  relation  specifying  the  coupled 
beam-circuit  system  in  the  backward  mode.  But  this  means  that  we 
have  only  three  independent  equations  to  determine  the  four  rf  variables. 
Under  these  circumstances  we  are  left  with  the  only  alternative  of  dis¬ 
regarding  one  of  the  rf  variables  in  (4-4),  say  the  beam  velocity  u  j, 
and  use  the  other  three  for  matching  at  the  pitch  change.  This  is  a  pro¬ 
cedure  which  has  been  used  in  the  past,  but  is  rather  questionable, 
because  the  forward  mode  is  bound  to  be  subject  to  some  reflections  at 
a  discontinuity  of  pitch  angle.  Only  in  the  very  special  case  that 
the  circuit  impedances  on  either  side  of  the  discontinuity  are  the  same. 
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will  there  be  no  reflections.  But  this  conditions  is  not  likely  to  occur 
in  the  general  case.  We  are  therefore  forced  to  adapt  six  independent 
boundary  conditions,  and  as  many  as  eight  conditions  if  pitch  reversals 
are  involved.  The  two  additional  conditions  in  the  latter  case  arise 
because  the  two  independent  m  =  -1  and  m  =  +1  beam  modes  are  both  present, 
traveling  along  the  beam  and  interacting  with  the  appropriate  space  har¬ 
monic  of  the  left-handed  and  right-handed  helix  sections,  respectively. 

From  Figure  4-5  it  is  noted  that  in  a  given  helix  section,  left-handed 
or  right-handed,  only  one  of  the  m  ■  -1  or  m  *  +1  beam-to-circuit  inter¬ 
actions  is  involved.  From  inspection  of  the  dispersion  diagrams  shown 
earlier  in  Figures  4-lb  and  4-2b  we  conclude  that  these  are  identical. 
Hence,  the  circuit  equations  of  the  m  =  -1  and  m  =  +1  interactions  will 
be  the  same.  However,  the  electronic  equations  for  m  =  -1  and  m  =  +1  are 
not  the  same  in  the  general  case.  Only  for  confined  focusing  are  they 
identical.  The  physical  reason  for  the  difference  is  the  presence  of 
the  dc  focusing  magnetic  field  which  introduces  a  nonsymmetric  effect 
with  regard  to  the  positive  (m  =  1)  and  negative  (m  =  -1)  angular 
dependence.  A  typical  example  is  the  Brillouin  focused  beam,  for  which 
the  propagation  factors  of  the  m  =  -1  and  m  =  +1  modes  are  different. 

Hence,  a  section  of  reversed  pitch  essentially  retains  the  same  circuit 
equation,  but  requires  a  difficult  electronic  equation,  except  for  con¬ 
fined  flow.  The  case  of  confined  flow  is  therefore  considerably  simpler 
to  analyze. 

4.4  DISPERSION  RELATIONS  FOR  THE  FORWARD  AND  BACKWARD  WAVES 


A  detailed  analysis  of  the  coupled  circuit  and  electronic  equations  leads 
to  the  same  form  of  the  dispersion  relation  for  the  forward  and  the 
backward  waves,  i.e.,  for  m  =  0  and  m  =  ±1.  The  relation  is  expressed 
in  the  matrix  form: 


0  = 


..U> 

J  "  '  \ 

z  1^1. 

cm  \  v  /  m 
'  cm 
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.  /  vo  v  vo 
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-B(i)f(i) 
m  m 

2  m  m 
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V‘>2 

p  m 

o 

0 

-1 

:(>  -  <■“’) 
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i  =  1,2,3 
m  =  0,  t 1 


(A. 5) 


The  parameters  appearing  in  the  equation  are  the  following: 

,(D 


B 


(i) 


,(i) 


- —  =  normalized  propagation  factors  of  the  three 

D 

e 

modes,  corresponding  to  i  =  1,2,3 
plasma  reduction  factors  for  the  same  modes 


f^^  =  coupling  coefficients  for  the  same  modes 

in 


v  =  phase  velocity  of  the  circuit  for  the  mth  space 
harmonic 


X  =  circuit  impedance  of  the  mth  harmonics 


cm 


m  R0 

1  ~  j  ( —  1 )  - - - —  =  loss  parameter 

cO^cO 
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where  is  a  diagonal  matrix  containing  the  three  propagation  factors  of 
the  mode  m. 

These  are  obtained  from  the  determinantal  equation  associated  with  (4.5) 

A  =  0  m  =  0,  il  (4.9) 

m 

Furthermore,  the  vector  A  in  (4.8)  contains  the  normal  mode  amplitudes 

__m 

of  the  three  component  modes,  which  are  independent  of  the  axial  coordi¬ 
nate  Zq  within  a  uniform  helix  section.  But  they  undergo  changes  at  any 
discontinuity  of  the  helix. 

ith  reference  to  the  variables  listed  after  (4.5)  one  should  note 
that : 

1.  The  plasma  reduction  factors  and  the  coupling  coefficients  are 
not  the  same  for  the  m  =  0  and  the  m  =  ii  modes. 

2.  The  circuit  impedance  Z  is  positive  for  forward  mode  m  =  0 
and  negative  for  the  backward  modes  m  *  il. 

3.  The  loss  parameters  x  is  the  same  for  m  =  0  and  m  =  il  except 

m  r 

for  the  sign  of  the  imaginary  part. 

With  these  restrictions  the  same  general  equations  apply  to  all  three 
modes  m  a  0,  and  m  **  il.  Note  that  in  this  formulation  all  the  modes 
have  positive  propagation  constants  ,  because  the  phase  velocities 

ITt 

are  all  positive.  The  circuit  part  of  the  m  =  il  modes  relate  to  the 
space-harmonic  voltage  and  current  v+1  and  i  ,  appearing  in  (4.8).  In 
the  next  step  one  must  convert  from  these  variables  to  the  corresponding 
fundamental  components  v^  and  i^  appearing  in  the  diagrams  of 
Figure  4-5. 
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4.5  SPECIFICATION  OF  THE  RF  VARIABLES  IN  THE  BACKWARD 
TRAVELING  WAVE 


In  view  of  the  strongly  coupled  space  harmonics  of  the  helix  this  pro¬ 
cedure  is  fairly  straightforward.  Let  us  first  look  at  the  propagation 
factor  of  the  m  =  0  backward  traveling  mode.  From  the  general  formula 
(4.3)  we  obtain 


B(i)  =  B(0 

0, backward  b 


B(i) 


(4.10) 


where  6  d  is  the  normalized  helix  pitch, 
e 


The  corresponding  backward  circuit  voltage  of  the  fundamental  component 
is  obtained  by  noting  that  the  total  rf  circuit  power  P  in  all  harmonics 
is  given  by  the  formula 


cm' 


2  P 


0,  il,  ±2, 


Hence 
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m 

2 

IB 

7 

cm 

< 

cr 

co 

(4.11) 


or 


m  =  fl  ,  (4.12) 


which  relates  the  normalized  rf  voltages  of  the  m  -  0  mode  and  the 
m  “  il  modes. 
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The  inverse  relation  applies  to  the  corresponding  rf  currents.  Since 


(4.13) 


we  obtain 

m  «  ±1  (4.14) 

In  this  backward  mode  Z  _  as  well  as  z  must  be  taken  negative,  but  we 

cu  cm 

have  used  absolute  signs  in  (4.12)  and  (4.14)  in  order  to  avoid  errors. 

Using  the  established  relations  (4.1 0) —(4.1 4)  between  the  m  *  0  and 
m  =  il  components  in  the  backward  mode,  the  relevant  rf  variables  for 
this  mode,  to  be  used  in  the  matching  procedure,  are  specified  by  the 
relation: 
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m  =  0,  ±1 

(4.16) 


By  the  described  procedure  we  have  been  able  to  express  the  rf  variables 
of,  the  backward  traveling  wave  in  the  form  (4.15)  which  involves  the 
m  =•  0  circuit  rf  variables.  This  is  a  necessary  step  in  the  matching 
process,  as  discussed  earlier. 


Before  concluding  this  section  on  the  backward  traveling  wave  we  must 
pay  proper  attention  to  the  discrimination  of  left-  and  right-handed 
helices.  In  light  of  the  earlier  discussions  it  follows  that  the 
m  *=  +1  and  m  =  -1  waves  do  not  simultaneously  exist  with  their  full 
range  of  composite  node  solutions,  described  by  the  superscript  i.  The 
full  set  corresponds  to  three  modes,  with  i  =  1,  i  =  2,  and  i  =  3, 
respectively.  But  either  the  m  =  +1  or  the  m  =  -1  wave  has  only  two 
composite  modes.  The  reason  for  this  is  that  one  of  the  two  m  =  ±1  set  of 
beam  waves,  both  of  which  have  about  the  same  phase  velocity,  is  way 
out  of  synchronism  with  the  corresponding  space  harmonic  component  of 
the  circuit.  We  can  easily  convince  ourselves  that  this  is  the  case 
by  studying  the  circuit  dispersion  diagrams  of  Figures  4-lb  and  4-2b. 
Hence,  for  either  m  =  +1  or  m  =  -1,  the  interaction  with  the  circuit 
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disappears,  leaving  us  with  the  two  decoupled  space  charge  modes  corre¬ 
sponding  to  the  particular  noninteracting  m  -  value. 

More  specifically,  for  a  left-handed  helix  the  m  =  +1  circuit  space  har¬ 
monic  is  way  out  of  synchronism  with  the  beam  (see  Figure  4- lb)  and, 
therefore,  the  m  =  +1  space  charge  waves  are  completely  decoupled  from 
the  circuit.  Hence,  in  the  left-handed  helix  section  only  two  modes 
corresponding  to  m  =  +1  exist.  Conversely,  in  a  right-handed  section 
there  are  only  two  modes  associated  with  the  m  -  -1  wave. 

These  considerations  imply  that  in  any  one  helix  section  we  have  a  total 
of  eight  independent  normal  mode  amplitudes,  which  will  be  discussed  in 
more  detail  in  the  next  section. 

4.6  SUMMARY  OF  NORMAL  MODE  AMPLITUDES 


In  general,  the  normal-mode  amplitudes  refer  to  the  independent  and  arbi¬ 
trary  components  in  the  normal  mode  vector  Am>  appearing  in  (4.8)  and 
(4.15).  In  this  form  they  are  arbitrary,  because  the  normal  mode  sc  u~ 
tions  have  not  yet  been  matched  at  the  discontinuities  between  the  helix 
sections.  In  effect,  the  matching  procedure  provides  exactly  sufficient 
conditions  to  determine  the  normal  mode  amplitudes  everywhere. 

Let  us  specify  explicitly  the  independent  normal  modes  in  left-handed 
and  right-handed  helix  sections,  respectively. 

4,6.1  Left-Handed  Helix 

The  fundamental  forward  wave  has  three  independent  modes,  corresponding 
to  m  =  0  and  i* =  1,  2,  and  3.  The  backward  wave  has  three  independent 
solutions  corresponding  to  m  =  -1,  and  two  solutions  corresponding  to 
m  =  +1 .  The  eight  independent  mode  amplitudes  can  be  summarized  as 
follows : 


23 


m  ■  0: 

A(1) 

0 

a(2) 

Ao 

a(3) 

Ao 

(4.17) 

m  ■  -1 : 

A(1) 

-1 

a(2) 

-1 

a(3) 

-1 

(4.18) 

m  “  +1 : 

A0) 

A+1 

a(2) 

A+1 

(4.19) 

4.6.2  Right-Handed 

Helix 

In  this  case  the  co 

t 

rresponding  normal  mode 

amplitudes  are: 

m  =  0 : 

A0) 

A0 

a(2) 

Ao 

a(3) 

Ao 

(4.20) 

*  m  =  -1 : 

A(1) 

-1 

a(2) 

-1 

(4.21) 

m  =  +1 : 

A(1) 

+  1 

a(2) 

A+1 

a(3) 

+1 

(4.22) 

4.7  COMPACT  FORMULATION  OF  THE 

FULL, SET  OF 

RF  EQUATIONS 

From  the  foregoing  discussion  it  is  quite  apparent  that  the  matching  of 
eight  normal-mode  amplitudes  at  the  intersections  between  different  helix 
sections  is  a  fairly  complex  mathematical  task.  It  can  best  be  accom¬ 
plished  by  etablishing  a  systematic  8-dimensional  matrix  formalism. 

For  convenience  and  clarity  of  the  exposition  let  us  introduce  sub¬ 
matrices  in  the  state  matrix  Q'  defined  by  (4.16),  (4.6),  and  (4.7). 
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|  2  (4.23) 

m  =  0,  ±1 

As  indicated  in  (4.23)  P  and  W  are  3x2  submatrices,  except  for  the 

~m  ~m 

decoupled  modes  corresponding  to  (4.19)  and  (4.21).  In  the  latter  case 

P  and  W  are  2x2  submatrices. 

~m  ~m 

We  can  now  proceed  to  use  (4.8)  and  (4.15)  for  the  various  waves  and 
modes  and  construct  matrix  relations  for  the  over-all  eight  rf  variables 
in  a  particular  section.  These  look  a  bit  different  for  the  left-handed 
and  right-handed  sections. 

4.7.1  Left-Handed  Helix  Section 


The  following  matrix  relation  applies: 
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(4.24) 
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The  left  column  specifies  all  the  eight  rf  variables  invjlved  in  the 
left-handed  helix  section.  The  first  8x8  matrix  on  the  right  is  a 
generalized  state  matrix  for  the  over-all  system.  The  second  8x8 
matrix  is  the  diagonal  exponential  matrix  containing  all  the  eight 
propagation  factors,  and  the  last  column  represents  the  eight  inde¬ 
pendent  normal  mode  amplitudes.  Equation  (A. 24)  is  an  extension  of 
(A. 8)  from  four  dimensions  to  eight  dimensions.  The  shaded  blocks  in 
the  equation  correspond  to  the  noninteracting  m  =  +1  space  charge 
waves . 

f 

A. 7. 2  Right-Handed  Helix  Section 

The  corresponding  relation  for  the  right-handed  section  is  the  following: 
3  2  3  3  2  3 


(A. 25) 
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Equations  (4.24)  and  (A. 25)  form  the  very  basis  of  the  analysis  of 
backward  wave  instabilities  in  multisection  TWTs.  All  the  elements  in 
the  8x8  matrices  are  known  from  the  specified  data  and  evaluation  of 
the  previously  stated  equations. 

It  has  been  common  practice  in  the  past  to  analyze  backward  wave 
instabilities  using  various  simplifications  in  the  model  and  in  the 
mathematical  procedure.  Most  of  these  approximations  are  justified  only 
by  the  ensuing  mathematical  simplifications  and  tractability ,  and  are 
not  justified  from  physical  reasoning. *  The  more  common  approximations 
can  be  easily  obtained  from  the  general  formulae  above.  Let  us  have  a 
look  at  these. 

A. 8  SPECIAL  CASES  AND  APPROXIMATIONS 


t  If  there  are  no  pitch  reversals  along  the  entire  helix  structure,  only 
one  of  the  two  m  *  ±1  space  harmonics  is  involved  in  the  interaction. 

Let  us  assume  that  the  helix  is  left-handed,  in  which  case  (A.2A)  applies. 
In  this  case  the  m  =■  +1  wave  can  not  be  excited  in  any  section  and  is 
therefore  nonexistent.  Hence,  the  corresponding  normal  mode  amplitudes 
specified  by  A^j  are  zero  everywhere. 

A,,  H  0  (4.26) 


These  two  modes  are  the  decoupled  space  charge  modes  which  had  to  be 
excited  in  some  other  right-handed  section.  But  since  these  do  not  exist, 
the  modes  are  nonexistent.  In  (A. 24)  these  are  represented  by  the  last 
two  equations,  corresponding  to  the  shaded  areas.  By  removing  all  the 
shaded  blocks,  the  resulting  equation  becomes  of  sixth  order,  and 
represent  an  exclusive  lof t-handed  helix  TOT. 


Similarly,  if  all  sections  are  right-handed,  the  m  *  -1  waves  are  not 
excited,  in  which  case: 


(4.27) 


The  relevant  equation  is  now  (4.25)  with  all  the  shaded  blocks  removed. 


Inspection  of  (4.24)  and  (4.25)  after  removal  of  shaded  blocks  show  that 
both  take  exactly  the  same  form:  ' 


v0(z0> 


w 


S+l<z0> 


_ L 

forward 

_ 1 _ 

backward 

wave 

wave 

column 

column 

(4.27) 


Even  if  the  equations  are  the.  same  form  for  left-handed  and  right-handed 
helices,  the  elements  in  the  matrices  are  not  the  same  for  in  =  -1  and 
m  =  +1.  They  are  the  same  only  for  confined  flow,  but  not  for  general 
balanced  beams,  including  Bril louin-beams .  These  more  subtle  points 
were  discussed  earlier. 
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Backward  wave  instabilities  are  often  analyzed  disregarding  the  forward 
traveling  wave.  This  is  legitimate  if  there  are  no  discontinuities  in 
the  form  of  pitch  changes  along  the  helix.  However,  the  use  of  this 
approximation  for  analyzing  the  effect  of  discontinuous  changes  in 
pitch  angle  is  very  questionable,  and  at  best  gives  a  crude  approxima¬ 
tion  to  the  starting  condition. 

However,  for  the  sake  of  completeness,  let  us  discuss  how  this  approxi¬ 
mation  follows  from  (4.27)  by  putting  the  forward  traveling  wave  ampli¬ 
tude  Aq  equal  to  zero. 

AQ  =  0  (4.28) 


Hereby,  (4.27)  is  simplified  to: 


—  — 

—  — 

w 

1o(zo) 

-p*l 

uil(zo> 

“•1 

Stl <zo) 

L 

which  is  just  the  same  as  the  original  equation 


(4.15). 


(4.29) 


In  the  present  work  we  shall  not  use  any  of  the  approximations  dis¬ 
cussed  here,  but  retain  generality  by  applying  (4.24)  and  (4.25). 

4.9  THE  TRANSMISSION'  MATRIX  OF  A  UNIFORM  HELIX  SECTION 


The  introduction  of  transmission  matrices  for  the  sections  is  particu¬ 
larly  convenient  for  the  purpose  we  have  in  mind. 


First,  let  us  write  (4.24)  and  (4.25)  in  the  following  compact  forms: 


-->?eZ0 

b(z0)  *  Me  e  Ar 


(4.30) 


~j5rzo 

MZq)  =  Mr  e  Ar 


(4.31) 


where  the  meanings  of  the  terms  follow  from  comparison  with  (4.24)  and 
(4.25).  Since  both  equations  are  of  the  same  form,  we  can  drop  the 
subscripts  C  and  r  and  simply  write: 


hJzQ) 


=  M  e 


~J9zr 


(4.32) 


The  proper  subscript  can  be  added  once  the  pitch  direction  has  been 
spec i f ied . 

Let  us  consider  the  uniform  helix  section  shown  in  Figure  4-6  and  apply 
(4.32)  at  both  ends  of  the  section.  Recalling  that  the  normal  mode 
vector  A  is  constant  and  independent  of  z^  within  the  uniform  section, 
we  find,  by  elimination  of  A  from  the  two  equations: 

b(zQ2)  *  M  e"jPL  M-1  b(zQ1)  ,  (4.33) 


where 


L 


z02  '  Z01 


(4.34) 


The  equation  is  in  t ransm i s s ion  form,  specifying  all  the  rf  variables 
b(Zp1'l  at  the  output  end  by  the  rf  variables  b(z^j)  at  the  input  end  of 
the  section.  The  transmission  matrix  T  is  given  by 


T 


(A. 35) 


-  „  ^  M-> 


Hence,  (A. 33)  takes  the  simple  form 


£(z02>  - 


(4.36) 


We  have  also  use  for  the  inverse  relation 

”  r's'v  ■  (4-37> 


where 

I  =  T_1  (A. 38) 

(i ' 

is  the  inverse  transmission  matrix. 

A. 10  THE  OVER-ALL  TRANSMISSION  MATRIX  OF  THE  MULTISECTION  TWT 

By  the  transmission  matrix  formulation  (A. 35)  we  have  established  the 
necessary  tools  for  a  conceptually  simple  and  mathematically  elegant 
description  of  a  multisection  TWT.  The  configuration  is  shown 
schematically  in  Figure  A-7.  For  each  of  the  N  sections  the  trans¬ 
mission  relation  (A. 36)  applies.  Some  of  the  sections  may  be  left- 
handed  and  some  right-handed,  and  care  must  be  taken  to  apply  the  cor¬ 
responding  version  of  the  transmission  matrix,  evaluated  either  from 
(A. 30)  or  (A. 31),  in  conjunction  with  (A. 35). 


As  noted  earlier,  all  the  rf  variables,  represented  by  the  columns  bCz^) 
and  bCz^)  in  (A.  36)  are  continuous  across  the  intersections.  The 
transmission  properties,  expressed  by  this  equation,  can  therefore 
easily  be  extended  to  comprise  any  number  of  sections,  and  in  particular, 
the  full  length  of  the  tube.  The  following  result  is  obtained: 
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Figure  4-7  Schematic  configuration  of  a  multisection  TWT 
consisting  of  a  total  of  n  uniform,  but 
different  sections. 
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(4.39) 


'  '  '  Vl  b<*ll>  • 


where  b(z  ,,)  refers  to  the  output  end  of  the  last  section,  and  b(z,,)  to 
~  r.z  —  1 1 

the  input  end  of  the  first  section.  If  we  write  (4.39)  in  the  form 


~  nz 


'll' 


f4. 40) 


the  over-all  transmission  matrix  T  is  given  by  the  product  of  all  the 
section  transmission  matrices 


T  *  Vn-l  •  *  *  ?2T1 

The  transmission  relation  expressed  in  (4.40)  *s  a  complete  and  detailed 
description  of  the  TWT  characteristics,  including  forward  gain  and 
backward  wave  properties.  All  the  elements  of  the  transmission  matrix 
T  can  be  specified  or  evaluated  from  the  detailed  relations  in  the 
earlier  part  of  this  report.  In  the  next  section  we  shall  show  that 
only  four  of  these  are  necessary  for  describing  two-port  terminal 
relat ions . 

4.11  TWO-PORT  TERMINAL  DESCRIPTION  OF  THE  MULTISECTION  TWT 

Let  us  express  the  transmission  relation  (4.40)  in  the  following  more 
explicit  and  detailed  form: 
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S,l<20> 


(4.44) 


Equation  (4.42)  can  then  be  subdivided  into  two  it.ponent  equations: 


b  (z  .)  -  T..  b 

-c  n2  -11  ~c 

(zn)  + 

(4.45) 

b  (z  .,)  -  T  b 

e  n2  —21  ~c 

(2ll)  +  ~22  ;(zll) 

(4.46) 

But 

at  the 

input  end  all  the  rf  variables  in  the  beam 

are  identically 

zero,  because  the  beam  enters  the  helix  unmodulated. 

responding  vector  b  (z,.)  is  identically  zero. 

-ell 

Hence,  the  cor- 

be(zll) 

;  0  , 

(4.47) 

and 

(4.45) 

and  (4.46)  simplify  to: 

h,(:  ,)  “ 

— c  n2 

In  -b<<*n) 

(4.48) 

b  (z  ,)  - 

e  n2 

T21  bc(zn> 

(4.49) 

As  a  consequence  of  the  fact  that  the  rf  beam  modulations  are  zero  at  the 

input  end  of  the  tube,  we  have  been  able  to  obtain  nice  separate  trans- 

mission  matrices  for  the  circuit  variables  b  and  the  beam  variables  b  . 
-  — c  ~e 

Of  these,  the  circuit  relation  (4.48)  is  by  far  of  most  significance,  and 
we  shall  not  use  (4.49)  any  further. 

It  is  quite  remarkable  that  a  simple  relation  like  (4.48)  describes  all 
the  small  signal  terminal  properties  of  the  multisection  TWT,  including 
forward  gain  characteristics  and  backward  wave  instabilities.  All  these 
properties  are  contained  in  the  2x2  transmission  matrix  Tp,  which  is 
known  from  the  computer  calculations  of  the  over-all  transmission 
matrix  T. 

Hence,  the  discussions  of  TWT  properties  are  reduced  to  a  discussion  of 
the  properties  of  the  submatrix  Tp  .  Only  in  the  simplest  cases,  and 
certainly  not  for  multisection  TWTs,  can  we  obtain  exact  or  approximate 
analytical  expressions  for  the  elements  of  ^ .  But  for  the  computer, 
determination  of  Tp  is  a  well  defined  problem  and  a  straightforward 
procedure  using  the  general  formula  (4.41). 
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5.0  BACKWARD  WAVE  INSTABILITIES 


From  this  base  on  we  can  proceed  to  determine  the  conditions  for  backward 
wave  instabilities  which  are  contained  in  the  properties  of  the  T^ 
matrix  in  (4.48).  Let  us  express  the  matrix  explicitly  by  its  elements: 


•'ll  X12" 

Y  V 

'21  '22 


(5.1) 


and  remember  the  relevant  equation  (4.48),  repeated  here  for 
convenience : 


W  ‘  ^c(zU>  (5'2> 

5.1  THE  INHERENT  BACKWARD  WAVE  INSTABILITY 

The  usual  argument  for  determining  backward  wave  instabilities  is  the 

following:  If  the  rf  circuit  variables  b^z^)  are  zero  at  the  output 

end,  but  nonzero  at  the  input  end,  i.e.,  bc(z^)  4  0,  then  the  backward 

gain  is  infinite  and  the  solution  unstable,  giving  rise  to  oscillations 

at  some  frequency.  This  condition  is  easily  obtained  from  (5.2)  by 

putting  b  (z  „)  SO,  which  results  in  the  equation: 

~c  n/ 

X  bc(zn)  =  0  (5.3) 


The  solution  of  this  homogeneous  set  of  equations  specifies  the  condi¬ 
tion  for  backward  wave  instabilities.  In  order  for  b  (z,,)  to  be  nonzero, 

~c  II 

the  determinant  of  (5.3)  must  be  zero 


X'  “  'll '22  12 '21 


0  , 


which  is  the  instability  criterion. 


(5.4) 
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The  elements  in  the  V -matrix  are  functions  of  all  the  relevant  param¬ 
eters  of  the  multisection  TWT.  It  reflects  the  properties  of  all  the 
sections  as  far  as  fixed  design  parameters  are  concerned,  such  as  helix 
impedances,  pitch  angles,  etc.  Moreover,  \  is  also  a  function  of 
operational  parameters,  in  particular  the  dc  beam  current  and  the  fre¬ 
quency  w.  We  are  not  at  liberty  to  specify  the  oscillation  frequency 
u,  which  is  an  unknown  parameter  together  with  the  starting  current  lfl  . 

Since  all  the  elements  in  the  \ -matrix  are  complete  quantities,  the 
determinantal  equation  (5.4)  contains  two  independent  equations,  repre¬ 
senting  the  real  and  imaginary  parts: 

|  .(u:,I0)  |  =  0  (5.5) 

r 

!‘(w,In)|  -  0  (5.6) 

i 

These  two  equations  define  implicitly  two  functional  relationships 
between  and  u>,  which  have  to  be  satisfied  simultaneously.  In 
Figure  5-1  these  two  relations  are  represented  by  the  intersecting 
curves  in  the  ui-Iq  plane.  The  Intersection  is  specified  by  point  A  in 
the  diagram  and  represents  the  starting  current  Iq^  and  oscillation  fre¬ 
quency  for  backward  wave  instabilities. 

The  numerical  procedure  would  involve  appropriate  search  routines  for 
determining  the  intersection  point  A  between  the  two  curves. 

At  this  stage  it  is  appropriate  to  point  out  some  common  errors  made 
in  determining  the  starting  current.  It  is  not  unusual  that  the 
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starting  condition  is  established  by  considering  only  one  equation, 
say  (5.5),  or  some  linear  combination  of  (5.5)  and  (5.6).  Regardless 
of  details,  the  point  is  that  only  one  relation  between  Iq  and  io  is 
specified.  In  Figure  5-1  this  would  correspond  to  disregarding  one  of 
the  two  curves,  say  the  one  corresponding  to  |x(cdq,Iq)|^  =  0.  In  this 
erroneous  procedure  one  is  then  plotting  the  remaining  curve,  which  in 
this  example  would  be  |'X(to,I-)|  =  0,  and  specifying  the  starting  con- 

dition  as  the  minimum  of  the  curve.  In  this  example  this  would  be 
point  B,  corresponding  to  an  alleged,  but  incorrect  starting  current 
Iq  and  frequency  u/. 

But  even  if  we  avoid  this  incorrect  procedure,  the  described  method 

using  both  equations  (5.5)  and  (5.6)  is  not  the  best  way  of  approaching 

the  instability  problem.  It  is  sure  that  the  method  gives  the  correct 

starting  condition  for  the  inherent  backward  wave  instability,  which 
*!  — — — 

has  to  do  with  internal  feed-back  loops  in  the  circuit-beam  system. 
However,  the  method  is  unable  to  cope  with  more  general  configurations 
characterized  by  additional  feedback  from  reflections  at  the  input  and 
output  terminals.  Expressed  differently,  the  described  procedure  is 
valid  only  for  perfectly  matched  input  and  output  terminals.  The 
reason  for  this  deficiency  is  that  the  relation  (5.2)  is  not  yet 
expressed  in  terms  of  its  forward  and  backward  traveling  rf  components 
and  the  corresponding  forward  and  backward  gain.  In  the  following  we 
shall  develop  a  procedure  along  these  lines,  and  arrive  at  a  better 
method  for  determination  of  backward  wave  instabilities. 

5.2  RESOLUTION  IN  FORWARD  AND  BACKWARD  COMPONENTS  IN  THE  OVER- ALL 
TRANSMISSION  SYSTEM  " 

As  already  pointed  out,  the  x_matrix  in  (5.2)  contains  all  the  informa¬ 
tion  on  the  terminal  behavior  of  the  multisection  TOT.  An  equivalent  way 
of  expressing  this  fact  is  the  statement  that  the  behavior  is  described 
by  the  eigenvalues  and  eigenvec tors  of  x • 


AO 
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The  elements  in  the  vectors  y^  and  y^  must  satisfy  certain  sign  require¬ 
ments  which  determine  the  choice  of  sign  in  (5.9).  The  details  are  dis¬ 
cussed  at  the  end  of  this  section. 

The  mathematical  procedure  is  now  to  express  the  circuit  rf  variables, 
which  are  given  by  the  vectors  b  .(z^)  and  bc(z^)  m  (5.2),  as  super¬ 
positions  of  the  two  eigenvectors  r^  and  r^.  Let 

5c(zn,)  =  Cf(zn2)  rf  +  Cb(zn2)  rb  (5.13) 
bcUll)  =  Cf(2ll)  rf  +  Cb<2ll)  rb  (5'l4) 


C,  (z  ,) 
b  n  2 


Ybcbczn) 


(5.17) 


What  is  the  physical  significance  of  the  mathematical  procedure  of 

expansions  into  eigensolut ions?  It  describes  the  resolution  of  the 

over-all  rf  circuit  variables  v^  and  i^  into  its  forward  and  backward 

traveling  components.  These  are  described  by  their  amplitudes  C,  and  C. 

f  b 

which,  of  course,  are  different  at  the  input  and  output  terminals. 


43 


Apparently,  the  gain  of  the  forward  and  backward  traveling  waves  are 
given  by: 


(voltage 

gain)  =  y£ 

(5.18) 

(.voltage 

(5.19) 

Hence,  the  respective  gains  are  specified  by  the  eigenvalue  vt  and  the 
inverse  eigenvalue  1 /y^ •  The  two  are  related  by  (5.10)  which  can  be 
expressed  as: 


(voltage  gain) 

(voltage  gain)  -  - ] — ; — r - -  (5.10) 

°  i  \  I 

This  equation  confirms  the  earlier  condition  (5.4)  for  backward  wave 
oscillations,  namely  |  \  j  *  0.  But  (5.20)  tells  us  more  than  (5.4).  If 
the  forward  gain,  at  the  oscillation  frequency,  is  zero,  or  very  small, 
the  procedure  of  putting  |  \  |  ■  0  is  bound  to  be  quite  inaccurate, 
because  the  right-hand  side  of  (5.20)  is  essentially  a  zero  over  zero 
expression. 

The  better,  and  correct,  procedure  is  to  determine  the  condition  for 
infinite  backward  gain  directly  from  (5.19),  which  is 

Yb  -  0  ,  (5.21) 


where  Is  specified  by  the  appropriate  expression  in  (5.9).  it 
follows  immediately  that  this  equation  gives 
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Hence,  it  looks  like  (5.21)  and  (5.4)  are  entirely  equivalent,  but  this 
is  only  from  a  superficial  viewpoint.  In  any  numerical  iterative  search 
routine,  the  end  result  is  a  small,  but  nevertheless  finite  value  of  y^ 
or  ] x  1  -  The  search  procedure  depends  on  the  variations  of  these  func¬ 
tions  around  their  zeros,  which  are  quite  different  for  y^  and  |x|  and 
favors  the  use  of  v,  .  The  discussion  in  conjunction  with  (5.20)  already 

'  D 

emphasized  this  point.  Moreover,  the  small,  but  finite  value  of  y^ 
obtained  as  the  end  result  of  the  search  procedure  also  tells  us  what  the 
backward  gain  is  under  these  circumstances,  namely  l/y^* 

The  forward  and  backward  solutions  are  basically  identified  by  the  nature 
of  the  eigenvectors ,  rather  than  the  eigenvalues.  The  latter  can  vary 
over  ;*  wide  range  because  the  forward  and  backward  gain  can  be  larger  or 
smaller  than  unity,  depending  on  the  frequency  and  other  operating 
parameters.  Hence,  under  genera.,  conditions  it  is  not  easy  to  identify 
the  forward  and  backward  components  from  the  two  eigenvalues . 

The  forward  wave  is  characterized  by  positive  power  flow  which  requires 
a  positive  real  part  of  the  impedance.  In  the  eigenvector  this  require¬ 
ment  is  reflected  in  the  signs  of  the  two  vector  elements.  The  real 
part  of  the  elements  must  have  the  same  sign. 

Conversely,  the  eigenvector  for  the  backward  wave  is  characterized  by 
opposite  signs  of  the  real  parts  of  the  vector  elements. 

The  eigenvectors  also  provide  us  with  direct  information  on  the  charac¬ 
teristic  Impedances  for  the  forward  and  the  backward  waves.  Thesa 
questions  are  discussed  in  Section  5.4.1. 

5 . 3  PROPER  PROCEDURE  FOR  DETERMINATION  OF  THE  INHERENT  BACKWARD 
WAVE  INSTABILITY" 

As  noted,  the  procedure  described  in  5.1  is  not  entirely  satisfactory. 

The  proper  procedure  is  based  on  the  condition  (5.2!)  rather  than 


(5.4).  The  real  and  imaginary  parts  are  both  zero.  For  convenience 
disregarding  the  subscript  b  we  then  have  the  two  equations: 

Yr(l0.  -r  +  J-t)  =  0  (5.23) 

Yi(1o'  “r  +  j“i>  "  0  *  (5.24) 

which  must  be  satisfied  simultaneously.  In  order  to  provide  a  better 
understanding  of  backward  wave  oscillations,  the  frequency  ..  is  assumed 
to  be  complex 


The  two  functional  relationships  (5.23)  and  (5.24)  are  sketched  in 
Figure  5-2  as  two  sets  of  intersecting  curves,  with  ^  constant  for  each 
curve.  The  intersections  between  curves  of  the  same  complex  frequency 
u,f  +  j-i  represent  the  solution  of  the  set  (5.23)  and  (5.24),  i.e., 
the  condition  for  inf inite  backward  wave  gain.  In  the  figure  the  solu¬ 
tion  is  given  by  the  dotted  line  B-C,  which  is  naturally  separated  in 
two  by  the  point  A  corresponding  to  a  real  frequency  Below  A  the 

solutions  have  posit ive  imaginary  frequency.  Apparently  this  corre¬ 
sponds  to  an  exponentially  decaying  solution  in  t lme ,  which  therefore 
is  stable .  The  solutions  above  the  point  A  have  negat lve  imaginary 
frequency  and  grow  uxpoent lally  with  time.  This  is  the  unstable  region. 
Hence,  if  we  visualize  that  the  current  Iq  is  increased  from  a  low 
value,  we  are  moving  upward  along  the  line  B-C.  The  system  is  nice  and 
stable  until  point  A  is  reached.  Point  A  repiesents  the  start ing  condi¬ 
tion,  which  is  characterized  by  a  rea I  oscillation  frequency  ,  ar.d  the 
starting  current  Iga-  If  the  current  is  increased  upwarus  t rum  A 
towards  C,  the  system  becomes  unstable,  with  exponentially  increasing 
rf  amplitude.  The  amplitude  is,  of  course,  limited  to  a  finite  value 
through  nonlinear  effects. 
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Figure  5-2  Sket'ch  illustrating  the  nature  of  backward  wave 

instabilities.  The  region  above  the  siiaded  outline  is 
unstable . 


5.4  THE  EFFECT  OF  MISMATCHED  INPUT  AND  OUTPUT  TERMINALS  ON  BACKWARD 
WAVE  INSTABILITIES 


In  the  preceding  treatment  of  the  backward  wave  instability  condition, 
it  was  assumed  that  the  input  and  output  terminals  were  matched.  This 
can  be  considered  to  be  the  basic  instability  criterion.  With  no  reflec¬ 
tions  at  the  input  and  output  terminals,  the  instability  is  due  entirely 
to  the  inherent  feedback  loops  in  the  backward  mode,  and  to  the  reflec¬ 
tions  taking  place  at  the  intersections  between  the  uniform  helix 
sections . 

In  this  section  we  shall  extend  the  analysis  to  the  more  general  con¬ 
figuration  of  nonmatched  input  and  output  terminals.  A  schematic  illus¬ 
tration  of  the  configuration  is  shown  in  Figure  5-3. 

In  a  practical  situation  one  would  like  to  determine  the  instability 
conditions  for  each  of  the  two  regions  which  are  separated  by  the 
attenuate-.  The  actual  configuration  is  shown  schematically  in  Figure 
5-4.  The  treatment  is  general  enough  to  be  applicable  for  this 
situat ion . 

In  the  previously  treated  inherent  instability  condition,  corresponding 
to  matched  input  and  output  terminals,  it  was  sufficient  to  consider 
equation  (5.1?)  describing  the  backward  wave  gain.  The  equation  (5.16) 
could  be  disregarded  because  of  the  forward  wave  amplitude  is  zero 
under  matched  conditions.  in  the  general  case  to  be  treated  here,  is 
not  zero.  It  will  be  convenient  to  introduce  the  ratio  C,/C,  at  the 

■  I  D 

input  and  output  as  new  variables.  But  first,  we  must  define  the 
char.ictet  1st  1c  impedances  of  the  forward  and  backward  waves. 

5.4.1  Forward  and  Backward  Characteristic  Impedances 

The  concept  oi  character  1st ic  impedance  iollows  directly  from  (5. 1 2)-(5. 14) . 
The  impedance  specifies  the  ratio  of  the  normalized  circuit  voltage  v 

o 


48 


ATTENUATOR 


Figure  5-3 
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Figure  5-4  Sketch  showing  the  TWT  separated  into  two  regions,  with 
the  attenuator  load  Z^t  serving  as  output  load  for  the 
first  region  and  Input  load  for  the  second  region. 


and  the  normalized  circuit  current  i^  for  each  of  the  two  modes  in  the 
absence  of  the  other.  The  normalized  impedances  follow  in  a  straight¬ 
forward  fashion  from  (5.11)  and  (5.12). 


The  forward  wave  characteristic  impedance  y^  is  given  by: 


(5.26) 


Note  that  j ^  is  a  normalized  impedance  because  Vq  and  i^  are  normalized 
with  respect  to  the  dc  voltage  and  current,  respectively. 

The  backward  wave  characteristic  impedance  f  is  given  by 


As  already  noted  in  Section  5.2,  the  real  parts  ot  f ^  and  are  posi¬ 
tive  and  negative,  respectively.  But  there  is  no  reason  to  expect  that 
yb  is  equal  to  minus  /^ .  Nor  are  these  impedances  specified  directly 
by  the  cold  circuit  impedances.  They  retlect  the  properties  of  the 
over-all  coupled  Deam  and  circuit  system  in  the  forward  and  backward 
directions,  respectively. 

5.4.2  The  Reflections  at  the  Input  and  Output  Terminals 

Introducing  the  ratio  of  the  forward  and  backward  mode  amplitudes  as  a 
variable,  we  can  use  (5.16)  and  (5.17)  to  obtain  the  following 
re  1  at  ion : 


n 


(5.28) 
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where  subscript  n  and  1  refer  to  the  output  and  input  terminals, 
respectively.  We  can  easily  determine  the  amplitude  ratios  from  (5.13) 
and  (5.14),  together  with  a  specification  of  input  and  output  load 
impedances  Z^j  and  Z  .  At  the  output  terminal 


where 


C 

o 


(5.29) 


(5.30) 


is  the  dc  beam  conductance. 

At  the  input  terminal  we  must  reverse  the  sign  because  the  rf  circuit 
current  is  defined  positive  in  the  positive  z-direction.  Hence 


Using  these  two 


( 

relat ions 

( 


together  with 


7. 


C 

o 


(5.13) 


nL  ('Q 

Vi¬ 


and 


(5.31) 


(5.14)  we  obtain: 


(5.32) 
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(5.33) 


It  is  clear  that  the  expressions  (5.32)  and  (5.33)  are  the  reflection 

coefficients  at  the  output  and  input  terminals,  respectively.  Denoting 

these  by  r>  and  .  we  have 
3  n  1 


'  n 


1 


n,  o 

h 


ZnL  Co 


A 


-  1 


(5.34) 


P 


1 


l  + 


Z1L  c,o 


A 


i  + 


Z1L  (,o 

h 


(5.35) 


One  should  note  that  the  reflection  coefficients  are  of  a  generalized 
nature  compared  to  reciprocal  transmission  systems,  because  the  charac¬ 
teristic  impedarces  are  different  in  the  two  directions. 

The  remaining  step  is  to  substitute  (5.32)  and  (5.33)  into  (5.28).  Also 
using  the  definitions  of  reflection  coefficients,  we  obtain  the  follow¬ 
ing  formula: 


V  n  ,‘b  '  1 


(5.36) 


This  is  the  new  instability  condition  which  is  characterized  by  unity 
loop  gain.  The  condition  of  unity  gain  is  clearer  i rom  inspections  o! 


Figure  5-5,  where  the  various  factors  in  the  product  (5.36)  are  repre¬ 
sented  in  the  closed  loop. 

Equation  (5.36)  can  be  expressed  in  the  equivalent  way: 

Yb  "  YfPlPn  =  0  (5-37) 

This  equation  represents  the  required  generalization  of  (5.21).  With 

matched  terminals,  r,  and  r  are  zero,  and  (5.37)  reduces  properly  to 

In 

(5.21). 

i 

If  we  define  a  modified  y^  by  the  relation 

Yb  *  Yb  “  Wn  •  (5‘38) 

it  follows  that  the  instability  condition  (5.37)  is  expressed  by 

-  0  (5.39) 

This  implies  that  the  general  nonmatched  configuration  can  be  treated 
by  exactly  the  same  mathematical  procedure  as  discussed  in  Section  5.3 
for  the  inherent  backward  wavt  instability.  We  simply  replace  y^  by  y£, 
defined  in  (5.38),  and  proceed  in  exactly  the  same  way. 

With  mismatched  input  and  output  terminals,  i.e.,  with  Oj  and  ,'n  both 
different  from  zero,  the  required  backward  gain  for  the  occurrence  of 
instabilities  is  less  than  infinity.  Hence,  the  starting  current  I  is 
expected  to  be  reduced  correspondingly. 

This  concludes  the  discussion  of  the  formal  structure  of  backward  wave 
instabilities.  We  have  established  a  general  mathematical  procedure  by 
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illustration  of  the  condition  (5.36)  of  unity 
as  criterion  for  backward  wave  oscillations. 
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which  very  general  configurations  of  multisection  TWTs  can  be  analyzed 

with  regard  to  both  forward  and  backward  wave  properties,  in  particular 

backward  wave  instabilities.  The  important  problems  of  optimum  design, 

i.e.,  the  design  of  multisection  TWTs  with  the  largest  possible  starting 

current  I„  without  undue  sacrifice  of  forward  wave  characteristics,  are 
us 

not  approached  analytically  in  the  present  work.  The  computer  program 
developed  from  the  theory  probably  represents  a  better  tool  than  any 
approximate  analytical  methods,  which  conceivably  could  be  developed  from 
the  theory. 
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'  6-°  PLASMA  reduction  factors  and  beam  coupling 

COEFFICIENTS  FOR  FORWARD  AND  BACKWARD  MODES 

The  dispersion  relation  and  state  vectors  for  the  coupled  beam-circuit 
system  in  a  uniform  helix  section  were  given  earlier  by  (4.5)— (4.7) . 

Among  the  parameters  appearing  in  these  relations  are  the  plasma  reduc¬ 
tion  factors  and  the  coupling  coefficients  f^\  where  m  =  0,  ±1 

m  m 

and  i  *  1 ,  2 ,  3 . 

A  considerable  amount  of  analytical  work  has  been  done  in  the  course 
of  this  program  to  determine  these  parameters  under  general  focusing 
conditions.  The  details  are  far  too  extensive  to  be  included  in  the 
final  report.  It  is  anticipated  that  a  full  account  of  beam  wave 
properties  under  general  focusing  conditions  will  eventually  appear  as 
a* separate  technical  report.  In  the  present  report  we  shall  state  a 
few  of  the  more  significant  results,  in  particular  the  plasma  reduc¬ 
tion  factors  and  the  beam-circuit  coupling  coefficients  for  confined  flow. 
We  shall  also  discuss  the  more  complex  Brillouin  focusing  condition  and 
the  corresponding  backward  wave  characteristics. 

6.1  BEAM  MODEL 

Plasma  reduction  factors  and  coupling  coefficients  are  concepts  that 
depend  on  certain  assumptions  concerning  the  distribution  of  rf  modu¬ 
lations  over  the  beam  cross  section.  The  details  of  the  longitudinal 
and  transverse  distributions  depend  in  large  measure  on  the  driving 
fields  from  the  surrounding  circuit.  However,  in  order  to  represent 
useful  concepts  in  TWT  design,  the  plasma  reduction  factors  and  coupling 
coefficients  should  be  largely  independent  of  the  details  of  the  circuit, 
except  for  a  minimum  of  geometrical  details  such  as  radial  dimensions 
of  the  beam  and  the  circuit. 
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This  implies  r.hat  one  has  to  make  reasonable  assumptions  concerning  the 
distribution  of  beam  modulations  over  the  cross  section.  In  the  more 
traditional  approaches  it  is  common  to  use  a  Bessel  function  distribu¬ 
tion,  which  is  simply  the  basic  term  in  an  infinite  normal  mode  expan¬ 
sion  of  space  charge  wave  components.  But  this  choice  is  dictated  more 
from  mathematical  convenience  rather  than  physical  reality.  In  view  of 
the  fact  that  the  edge  regions  of  the  beam  are  modulated  more  strongly 
than  the  central  part,  the  Bessel  function  distribution  is  obviously 
not  a  good  choice,  because  the  Bessel  function  has  its  maximum 

at  the  center  and  decays  towards  the  beam  edge,  i.e.,  the  variation  is 
exactly  opposite  from  the  actual  situation. 

With  this  in  mind  it  seems  logical  to  assume  the  very  simplest  distri¬ 
bution,  which  is  that  of  constant  velocity  and  displacement  distributions 
over  the  beam  cross  section.  Although  this  is  not  the  actual  physical 
dis  ribution  in  the  TWT,  it  is  a  better  choice  than  the  traditional 
Bessel  function  distribution. 

Accordingly,  the  longitudinal  and  transverse  dynamic  variables  for  the 
m  «  0  and  the  m  ■  £l  modes  are  specifed  by  the  schematic  diagrams  in 
Figure  6-1.  On  this  basis  we  have  developed  a  complete  field  theory 
for  determination  of  the  associated  fields  which,  in  turn,  specify  the 
plasma  reduction  factors  and  coupling  coefficients. 

6.2  CONFINED  FLOW 


Confined  flow  represents  a  limiting  case  of  focusing  which  is  never 
quite  achieved  in  a  practical  tube,  because  it  requires  infinite  mag¬ 
netic  field,  or  zero  beam  current.  But  the  condition  can  be  approached 
to  a  degree  which  justifies  the  use  of  this  concept. 

As  discussed  earlier,  confined  flow  is  by  far  the  simplest  from  a  con¬ 
ceptual  viewpoint.  Since  no  beam  rotation  takes  place,  the  m  ■  -1  and 
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6-1  Schematic  diagrams  of  the  assumed  distributions  of  the 
dynamic  variables  over  the  beam  cross  section. 

a)  Fundamental  mode,  ra  »  0 

b)  Backward  wave  modes,  m  ■  "1. 
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m  ■  +1  beam  modes  are  identical,  except,  of  course,  for  the  opposite 
angular  variations.  The  plasma  reduction  factors  and  the  coupling  coef' 
ficients  are  the  same  for  m  =  +1  and  m  =  -1.  But  these  differ  from  the 
corresponding  factors  of  the  fundamental  mode  m  ■  0. 


6.2.1  Plasma  Reduction  Factors  for  Confined  Flow 


For  the  fundamental  mode  m  =  0,  the  square  of  the  plasma  reduction 
factor  is  given  by  the  expression: 


o  I,(yO 

R2  „  -  1-2  A . K  G 

m=0 


OVF7  “01  <YV  Yrb) 


(6.1) 


The  variables  appearing  in  the  equation  are: 

Beam  radius,  r^ 

Helix  radius,  r  (6.2) 

3 

Propagation  factor,  y 

The  functions  I,(Yr,)  and  I_(Yr  )  are  modified  Bessel  functions,  and 
lb  U  a 

G_.  (yt  ,  yv.  )  is  defined  as  follows: 

U 1  3  D 

G01(Yra*  Yrb}  =  Kl(Yrb)  +  Ka(Yra}  1l(yTb)  (6*3) 

The  propagation  factor  y  basically  refers  to  the  propagation  factors 

of  the  three  forward  traveling  modes,  corresponding  to  i  =  1,  2,  and 
3.  But  these  are  sufficiently  close  to  justify  the  use  of  the  same  y 
for  all  three  components.  Hence,  for  practical  purposes 

Y  -  Be  (6. A) 
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,(yr  ,  yr.)  =  I,  (yr  J  Kn(yr.)  +  K.(yr  )  In(yr.) 


(6.6) 


The  plasma  reduction  factors  Rq  and  R+^  are  plotted  in  Figures  6-2 

and  6-3  for  several  values  of  the  beam-to-helix  diameter  ratio  r, /r  . 

b  a 

From  the  expressions  for  Rq  and  R+^  one  can  show  that 


lim 

yV0 


lim  R 

YV° 


0 


(6.7) 


lim  Rq  =  lim  R+^  =  1  (6.8) 

yr  -**>  yr,-*° 

b  b 


The  first  of  these  two  limiting  cases,  R  =  0,  corresponds  to  an 
infinitely  thin  beam  in  which  case  there  are  no  longitudinal  fields,  only 
transverse  field  components.  The  second  limiting  case,  R  =  1,  corre¬ 
sponds  to  infinite  normalized  beam  radius,  which  essentially  represents 
a  one-dimensional  case.  The  fields  are  pure  longitudinal,  and  the 
transverse  field  components  are  zero. 
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Figure  6-2  Plasma  frequency  reduction  factor  Rq  for  forward  wave 
interaction  and  confined  flow. 
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6.2.2  Coupling  Coefficients  for  Confined  Flow 


The  coupling  coefficient  f^  is  a  measure  of  the  strength  of  the 

tn 

coupling  between  the  beam  and  the  current.  The  coefficient  appears  in 
the  dispersion  relation  (4.5)  and  the  state  vectors  (4.6). 


For  the  fundamental  mode,  m  *  0, 
the  ratio  of  the  circuit  current 
the  beam 


Che  coupling  coefficient  is  defined  as 
1^  and  the  longitudinal  current  1^  in 


f 


0 


(6.9) 


The  coupling  coefficient  is  obtained  from  the  equation: 

£  .  21l(Yrb}  1 _ 

°  *  Yrb  10(Yra) 


(6.10) 


The  expression  is  plotted  in  Figure  6-4  for  several  values  of  the 

ratio  r,/r  . 
b  a 

The  mode  symmetry  in  the  backward  traveling  waves  is  such  that  the  cor¬ 
responding  coupling  coefficient  does  not  have  the  same  simple  physical 
interpretation  as  (6.9)  for  the  fundamental  mode  m  ■  0. 

With  reference  to  Figure  6- lb,  the  longitudinal  current  density  in  the 
m  -  il  modes  is  specified  by 


1  b ( r ,  •»,  i)  -  t^U)  e'-*'  p^r,  (6.11) 

where  i ^ ( z )  is  interpreted  as  the  one-dimensional  equivalent  of  the 
constant  current  density  in  the  fundamental  mode.  The  actual  current 
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density  is  proportional  to  the  radius  and  varies  azimuthally  as  exp  (ij6). 
The  total  integrated  current  is  therefore  zero.  The  entire  theory  in  the 
present  report  is  based  on  interpreting  ^(z)  and  (6.11)  as  an  equiva¬ 
lent  one-dimensional  current  density.  The  coupling  coefficient  f+  is 
defined  on  the  basis  of  this  variable,  in  the  following  way: 


f 


♦1 


i  (z)  r  do 


(6.12) 


where  i^(z)  is  the  surface  current  density  in  the  surrounding  circuit, 
in  the  m  »  ±1  mode. 

Evaluation  of  the  coupling  coefficient  f+^  yields  the  following 
expression: 


il 


1.(6  r) 
1  e  a 


10(eerb)  ‘ 


2Il(6erb> 

Vb 


(6.13) 


The  equation  is  plotted  in  Figure  6-5.  It  is  noted  that  f ^  can  be 
larger  than  unity.  This  is  simply  due  to  the  fact  that  the  coupling 
coefficient  for  the  antisymmetric  backward  wave  modes  does  not  have  the 
same  simple  physical  interpretation  as  the  corresponding  coefficient  f^ 
for  the  fundamental  mode.  In  particular,  the  numerical  value  of  f 
depends  on  how  we  define  the  "one-dimensional  equivalent  beam  current 
density"  i^(z).  The  definitions  used  in  the  present  work  are  stated  in 
(6.11)  and  (6.12),  and  all  remaining  equations  are  consistent  with  these 
definitions. 


6.3  BHILLOl'IN  FLOW 

In  the  present  program  an  extensive  analysis  has  been  made  oi  the 
nature  of  electron  beam  waves  in  a  general,  balanced  beam  model.  Tills 
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model  includes,  as  special  cases,  the  confined  flow  condition  treated 
in  the  previous  section  and  the  Brillouin  flow  condition.  The  dispersion 
relation,  and  thereby  the  plasma  reduction  factors  and  the  coupling 
coefficients,  are  specified  by  an  implicit  equation  of  sixth  order. 

For  the  simple  case  of  confined  flow  one  can  evaluate  these  factors 
explicitly.  However,  for  Brillouin  beams  the  situation  is  not  so  simple. 
The  roots  of  the  sixth  order  equation  must  be  solved  by  numerical  com¬ 
puter  procedures,  and  only  fairly  crude  approximations  can  be  obtained 
analytically.  , 

The  roots  of  the  sixth  order  equations  determine  a  set  of  six  beam 
waves.  In  addition  to  the  regular  space  charge  waves  the  beam  supports 
two  cyclotron  waves  and  two  quasisynchronous  displacement  waves. 

■■  The  analysis  has  been  carried  all  the  way  to  a  detailed  specification  of 
all  the  elements  in  the  sixth  order  dispersion  relation.  But  the  estab¬ 
lishment  of  numerical  solutions  and  interpretations  of  the  characteris¬ 
tics  of  the  beam  waves  would  represent  a  large  effort  much  beyond  the 
scope  of  the  present  program.  Therefore,  the  following  treatment  is 
limited  to  presenting  approximate  solutions  for  the  backward  wave 
characteristics  in  Brillouin  beams. 

6-3.1  Plasma  Reduction  Factors  for  the  Backward  Waves  (m  »  tl)  in 
Brillouin  Flow  ~ 


For  confined  flow  we  found  that  the  plasma  reduction  factors  for  the 
m  =  +1  and  m  =  -1  modes  were  the  same.  The  reason  for  this  result  is 
that  the  beam  motion  is  purely  longitudinal.  Hence,  for  reasons  of 
symmetry  there  can  be  no  difference  in  the  behavior  of  the  m  =  +1  and 
m  =  -l  modes,  except,  of  course,  the  opposite  azimuthal  variations. 

On  the  other  hand,  for  nonconfined  flow,  such  as  the  Brillouin  condi¬ 
tion,  the  m  =  +1  and  m  =  -l  modes  have  different  plasma  reduction 
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factors.  They  also  depend  on  the  direction  of  the  dc  magnetic 
focusing  field. 

An  approximate  perturbation  procedure  is  reasonably  accurate  for  the 
m  =  il  inodes  because  this  particular  symmetry  maintains  the  beam  cross 
section  everywhere  under  rf  conditions.  Hence,  the  transverse  dis¬ 
placement  is  not  likely  to  play  an  important  role  for  the  space  charge 
field.  On  the  other  hand,  for  Brillouin  flow,  the  circularly  sym¬ 
metric  mode,  m  »  0,  is  characterized  by  uniform  radial  expansions  and 
contractions,  which  clearly  contribute  to  the  axial  space  charge  field 
and  therefore  to  the  plasma  reduction  factor.  Accordingly,  we  limit 
the  approximate  procedure  to  the  m  »  il  backward  modes. 


The  results  for  the  m  -  +1  mode  can  be  stated  as  follows: 


Slow  wave: 


Fast  wave: 


+  1  ,s 


+  1  ,b 


il 


-  R  ,  i 


1 

S  2 
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il 


(6.14) 


(6. IS) 


where  R+1  is  the  plasma  reduction  factor  for  conf ined  flow,  specified 
by  (6.5)  or  Figure  6-3.  The  upper  signs  in  (6.14)  and  (6.15)  apply  for 
positive  direction  of  the  dc  magnetic  focusing  field,  the  lower  sign 
for  negat ive  direction  of  the  field. 

For  the  m  •  - 1  mode  the  con espondlng  equations  are: 


Slow  wave:  K  ,  ■  K  , 

-  1  ,s  i  l 


(6. 16) 


K-l.b  "  Kil  ‘ 


Fast  wave: 


(6. t  7) 


where  the  upper  and  lower  signs  refer  to  the  same  conditions  of 
focusing  as  stated  above. 

According  to  these  relations  the  plasma  reduction  factors  depend  on  the 
sign  of  m,  and  the  sign  of  > he  dc  magnetic  field.  The  following  general 
relations  apply: 


(6.18) 

Rf(m,BQ)  =  Rf(*m,  -BQ) 

(6.19) 

The  consequences  of  these  relations  for  backward  wave  oscillations  under 
Brillouin  flow  conditions  are  such  that  we  expect  the  starting  condition 
to  be  different  for  reversed  direction  of  magnetic  field.  Or  equiva¬ 
lently,  as  shown  by  the  general  relations  (6. 14)— (6. 19) »  a  right-handed 
and  a  left-handed  helix  do  not  have  the  same  starting  conditions  for 
the  same  direction  of  magnetic  field.  This  difference  is  expected  to 
affect  both  the  starting  current  and  the  frequency  of  oscillation, 
and  should  be  observable  experimentally  in  uniform  Brillouin  focused 
TWTs . 

In  a  periodically  focused  TWT ,  the  focusing  fields  in  alternate  sections 
correspond  to  the  positive  and  negative  Brillouin  field.  Accordingly, 
the  plasma  reduction  factors  alternate  between  the  upper  and  lower  values 
in  (6 . 14)- (6 . 1 7) .  The  synchronization  condition  for  backward  wave 
oscillations  also  changes  periodically  with  the  same  periodicity.  Due 
to  the  dominant  use  of  periodic  Brillouin  focusing  in  TWTs,  an  impor¬ 
tant  question  is  how  to  approach  the  problems  associated  with  the 
periodic  nature  of  the  electron  beam.  We  could  suggest  to  incorporate 
the  periodic  sections  into  the  over-all  forward  and  backward  inter¬ 
action  system  analyzed  in  Chapters  4  and  5.  In  principal,  one  could 
think  of  introducing  the  concept  of  a  "sectioned  beam"  in  the  same  way 
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as  done  for  sectioned  helices.  However,  this  is  more  difficult  for  the 
beam,  because  of  mode  conversions  taking  place  at  the  field  reversals, 
but  it  may  still  be  possible  to  implement.  It  would  require  more 
detailed  knowledge  of  all  the  six  beam  modes  discussed  earlier,  and  of 
the  appropriate  matching  conditions  of  the  field  reversals.  In  itself 
this  is  a  considerable  task  far  beyond  the  scope  of  the  present  program. 

In  lack  of  such  a  description  we  are  forced  to  accept  approximate 
procedures.  If  we  insist  on  using  one  single  plasma  reduction  factor 
for  the  entire  length  of  the  periodically  Brillouin  focused  beam,  the 
best  choice  appears  to  be  the  average  value  of  the  reduction  factors  in 
alternate  sections.  But  from  (6. 14)— (6. 17)  we  observe  that  the  average 
values  are  just  R+p  i.e.  ,  the  reduction  factor  for  confined  flow. 

Hence: 

0O  =  (R+  )  (6.18) 

Brill  ~  confined 

In  this  approximate  description  the  periodic  Brillouin  beam  looks  like 
a  confined  beam,  as  far  as  the  m  =  ±1  modes  are  concerned.  Physically, 
this  picture  is  quite  acceptable  in  view  of  the  fact  that  the  average 
dc  rotational  frequency  is  zero,  i.e.,  the  same  as  in  confined  flow. 

In  the  same  approximation,  we  would  expect  the  average  coupling 
coefficient  to  be  equal  to  the  coupling  coefficient  for  confined  flow: 

=  (f  +  i)  (6.19) 

Brill  “  confined 

This  concludes  the  discussion  of  backward  wave  characteristics  for  the 
Brillouin  focused  beam.  Further  analysis  is  required  of  the  special 
problems  of  periodic  Brillouin  focusing  and  its  implementation  in  the 
over-all  forward  and  backward  wave  system  described  in  Chapters  4  and  5. 


70 


LIST  OF  SYMBOLS 


a 


?b 


?f 


a 

b 


b 

„e 


c 


d 


Column  vector  specifying  the  rf  variables  along 
the  tube. 

Column  vector  for  the  backward  wave  mode. 

Column  vector  for  the  forward  wave  mode. 

Column  vector  for  the  mth  space  harmonic  component. 
Overall  column  for  the  rf  variables. 

f 

Circuit  part  of  b. 

Beam  part  of  b. 

Velocity  of  light  in  vacuum. 

Pitch  of  helix. 


et  Charge  of  the  electron. 

f^  Coupling  coefficient  of  the  ith  mode  of  the  mth 

m  space  harmonic. 


i^  Backward  wave  mode  circuit  current  normalized  to 

the  dc  current  (mode  m  =  0). 


Backward  wave  mode  circuit  current  normalized  to 
the  dc  current  (mode  m) . 


Forward  wave  mode  circuit  current  normalized  to  the 
dc  current  (mode  m) . 


m 


Space  harmonic  number. 


m  Rest  mass  of  electron, 

o 

m?  Relativistic  mass  of  electron  in  the  longitudinal 

direction. 


m. 


Relativistic  mass  of  electron  in  the  transverse 
direction. 


The  ith  component  of  the  state  vector  for  the  mth 
space  harmonic. 
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A 


~m 


Normal  mode  vector  of  the  mth  space  harmonic. 

Normalized  propagation  factor  of  the  ith  mode  of 
the  mth  space  harmonic. 


Normal  circuit  mode  amplitude  of  forward  mode  in 
the  TWT . 


C,  Normal  circuit  mode  amplitude  of  backward  mode  in 

L'  the  TWT . 


D 

_m 


Diagonal  matrix  containing  the  normalized  propaga¬ 
tion  f actors  . 
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LIST  OF  SYMBOLS  (CONTINUED) 


-*  ,-*s 
E  (r ) 
m 

Goi(8V6rb> 


W 

W 


mb  ) 

o  er 


W 


Electric  field  of  the  mth  space  harmonic  component. 
Combination  Bessel  function. 

DC  beam  current. 

Inverse  of  transmission  matrix  T. 

Modified  Bessel  function  of  first  kind  and  zero  order. 

Modified  Bessel  function  of  first  kind  and  first 
order. 

Modified  Bessel  function  of  second  kind  and  zero 
order. 

Modified  Bessel  function  of  second  kind  and  first 
order. 

Normalized  length  of  one  helix  section. 

Overall  state  matrix  for  left-handed  helix. 

Overall  state  matrix  for  right-handed  helix. 

Submatrix  of  state  matrix  O'. 

..m 

State  matrix  of  the  mth  space  harmonic  component. 

Fart  of  O' . 

_m 

Plasma  reduction  factor  of  the  ith  mode  of  the  mth 
space  harmonic  component. 

Normalized  displacement  of  the  mth  space  harmonic 
component . 

Overall  transmission  matrix. 


T(P) 


Transmission  matrix  of  section  p. 
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U 

m 


V 

o 


W 

_m 


z 

o 


cm 


n2 


'11 


'12 


nL 


o 

K 

m 


fcl 

c 

m 


Submatrices  in  the  overall  transmission  matrix  T. 


Normalized  rf  beam  velocity  of  the  mth  space 
harmonic  component. 

DC  beam  velocity. 

Submatrix  of  O’. 

~m 

Normalized  length. 

Circuit  impedance  of  the  mth  space  harmonic,  referred 

to  the  circuit  location  r  ■  r  . 

a 

Normalized  position  at  output  end  of  last  section. 
Normalized  position  at  input  end  of  first  section. 
Load  impedance  at  the  input. 

Load  impedance  at  the  output. 

Propagation  factor  of  the  mth  space  harmonic. 
Relativistic  factor. 

Eigenvalue  of  the  ^-matrix. 

Eigenvalue  of  the  x-matrix. 

Permittivity  of  free  space. 

Lo«<«  parameter  of  the  mth  space  harmonic  component. 

Reflection  i oof fi' lent  at  the  input  end  of  the  TLT. 

Reflation  coefficient  at  th'-  the  output  end  of 
the  TVT . 
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LIST  OF  SYMBOLS  (CONTINUED) 


Angular  frequency. 

Complex  angular  frequency. 

Imaginary  angular  frequency. 

Non-reiativistic  plasma  frequency. 

Relativistic  plasma  frequency. 

Plasma  frequency  normalized  to  the  operating  frequency. 
A  2  x  2  submatrix  equal  to  T^. 


The  elements  of  x* 

Determinant  of  x* 

Real  part  of  jx|  . 
Imaginary  part  of  |X|  • 
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